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Bottom and Charm quark masses from lattice NRQCD 

NRQCD collaboration presented by K. Hornbostel a * 

a Dept. of Physics, Southern Methodist University, Dallas, TX 75275 



We present new values for the MS masses of b and c quarks based on lattice NRQCD simulations of the T(bb) 
and ip(cc) systems. These include three measurements of the b mass based on quenched simulations with lattice 
spacings ranging from 0.05fm to 0.15fm, which we find to be largely independent of lattice spacing. In addition, 
we find a consistent value from an unquenched simulation at 0.08fm. 



1. Simulations 

The T(bb) and ip(cc) systems possess several 
properties which permit accurate lattice simula- 
tions. They are physically small, allowing the 
use of small volumes. They are insensitive to 
the presence of light quarks and to uncertain- 
ties in their masses. They are well understood 
phenomenologically, which aids in the estimation 
of systematic errors. Their spin-averaged mass 
splittings are insensitive to the b or c masses, al- 
lowing independent tuning of the lattice spacing 
and bare masses. Their properties have been pre- 
cisely measured. Finally, they may be efficiently 
simulated using a nonrelativistic effective action, 
NRQCD. 

To compute the spin-averaged spectra we em- 
ployed an NRQCD action for the heavy quarks, 
improved to include leading and next-to-leading 
order corrections for relativity and discretization 
errors. We used the average plaquette to tad- 
pole improve the link operators. Several collab- 
orations generously provided gauge held configu- 
rations for both zero and two light quark flavors, 
generated using the standard (unimproved) Wil- 
son action for gauge fields and staggered fermions 
for the light quarks. We removed perturbatively 
the leading discretization errors in the spectrum 
due to use of this unimproved gluonic action |Q . 
The details of our simulations and references for 
the gauge field configurations appear in Ref. H. 



*Work done in collaboration with C. T. H. Davics, 
G. P. Lepage, C. J. Morningstar, J. Shigemitsu and 
J. Sloan. Support provided by the DOE, NSF, PPARC, 
the Leverhulme Trust and the Fulbright Commission. 
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Figure 1. MS 6-quark mass M&(Mb) at various 
lattice spacings a. For rif = 0, (3 = 6.2, 6.0 and 
5.7, respectively. For nf = 2, /3 = 5.6. 



2. Extracting M b and M c 

Results from these simulations allowed us to 
obtain a value for the c quark MS mass M c and 
values for the b quark mass Mj at several lattice 
spacings, as follows. We first measured the lattice 
spacing a for each simulation by comparing either 
the computed spin-averaged IP-IS or 2S-1S split- 
tings with data. Because these splittings are very 
insensitive to the bare mass aM° , we were able to 
determine a without first finely tuning the mass. 
The results were accurate to within about 5% for 
the T system, 15% for tp. 

Using this a, we tuned the bare mass aM° so 
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Table 1 

Perturbative parameters connecting the bare mass M° to the MS mass M(n). 



a 


7lf 


aM 6 ° 


a\x 




aq* 


a r (q*) 


5.7 





3.15 


5 


-1.18 


3.03 


.1911 


6.0 




1.71 




-1.63 


5.73 


.1315 


6.2 




1.22 




-1.79 


9.80 


.1088 


5.6 


2 


1.80 




-1.61 


5.34 


.1570 






aM c ° 










5.7 





0.80 


7 


-2.85 


23.02 


.1086 



that the kinetic mass M^ ni determined by fitting 
to the dispersion relation 



E(p) = £7(0) 



P 



P 



2M k 



(1) 



agreed with experimental values for the T or r\ c . 

Once the fundamental parameters aM° and a 
were fixed, we were able to convert quantities to 
physical units. By using the expression 



M° = aM 



o ( M > 



cxpt 



\aM kil 



(2) 



to express the bare mass in GeV, we minimized 
sensitivity to uncertainties in aM and a. The 
bulk of aMki rl is made up of 2aM°, with dynamics 
providing the remainder. As a result, the ratio 
aM JaM^a is very insensitive to errors in the 
dynamics from uncertainties in aM° and a; these 
represent only an error in an order v 2 correction 
to 1/2. As an extreme illustration, computing 
Mj using a value for oMj as far off as aM® gives 
a result off by only 25%. In contrast, using a -1 
directly in the form M° — (aM°)a~ 1 rather than 
Eq. ^| would unnecessarily promote errors in a -1 
by one power of the mass. 

Refs. [| and WJ give perturbative expressions 
relating the MS mass M and the bare mass M° to 
the pole mass, respectively. We combined these 
to obtain 

. a P (q*) 



M(n) = M° 



1 



(3) 



with parameters listed in Table [l| We used the 
scheme in Ref. || to determine the scale q* for 



a P in this relation. In particular, we chose to 
extract M(/x) at a scale \i such that q* was well- 
defined. Values for a P came from measurements 
of the plaquette in these same simulations. 

Finally, we ran M(fi) to its own scale using 
three-loop evolution. Our preliminary results are 
presented in Table |^. Insensitivity of the final 
result to the bare mass aM® is evident. The 
systematic errors quoted account for higher-order 
relativistic and discretization corrections, most of 
which were based on estimates from a potential 
model. They do not account for errors associated 
with using the wrong number of light quarks. Val- 
ues for Mb at three lattice spacings and for zero 
and two light flavors rif also appear in Figure [p. 

We also applied an alternate method to deter- 
mine the T kinetic mass. Rather than fitting 
to the dispersion relation above, we determined 
aAfki n from 



aM] 



kin 



2 Zr„ aM" + aE 



NR 



2a£7 



(4) 



The first term on the right hand side gives twice 
the pole mass. The next two represent the energy 
for this state obtained in the simulation with the 
self energy of the quarks subtracted, which gives 
the binding energy. The quark mass renormal- 
ization Z rn and self energy aEg were computed 
perturbatively to one loop in Ref. . 

We used this kinetic mass to convert the bare 
mass to physical units and obtained the MS 
masses as above. Results were consistent with 
the first method, but with larger errors due to 
unknown higher-order terms in Z m and aEg. 
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Table 2 

MS masses for b and c quarks (preliminary). The lattice spacing a is determined from spin-averaged 
splittings between IP and IS states. Values for M° are from Eq. 0, rather than from (aM°) a' -1 ', as 
discussed in the text. Errors are statistical and estimates of higher order relativistic and discretization 
corrections, respectively. Systematic errors are combined in M° and M(M). The final error in M(M) 
estimates neglected higher-order perturbative contributions. We do not include an error estimate to 
account for quenching. 



a 


Tlf 


aM° 


a[-% (GeV) 


4;i\ s (GeV) 


M° (GeV) 


M b (M b ) (GeV) 


5.7 





3.15 


1.41(4)(2)(5) 


1.36(13) (2) (4) 


4.22(5)(8) 


4.15(5)(8)(3) 


6.0 




1.71 


2.59(5)(3)(1) 


2.45(8)(3)(1) 


4.11(3)(2) 


4.28(3)(3)(3) 






1.80 






4.16(3)(2) 








2.00 






4.21 (3)(2) 








3.00 






4.32(3)(2) 




6.2 




1.22 


3.52(14) (4) (0) 


3.24(15) (4) (0) 


3.99(6)(2) 


4.31(6)(3)(3) 


5.6 


2 


1.80 


2.44(6)(3)(1) 


2.38(10) (3) (1) 


4.08(4)(3) 


4.26(4)(3)(5) 






aM» 








M C (M C ) 


5.7 





0.80 


1.23(4)(3)(8) 


1.20(20) (3) (8) 


.98(3)(9) 


1.20(4)(11)(2) 


errors: 






stat, rel, discr 




stat, sys 


stat, sys, pert 



3. Conclusions 

We presented new results for M b and M c from 
lattice NRQCD calculations of the T and ip spin- 
averaged spectra. Our results for M b showed lit- 
tle dependence on a or on the presence of light 
quarks. They are consistent with those pre- 
sented at Lattice 97 and summarized in Ref. ||: 
M fc = 4.15(5)(20) and M c = 1.525(40) (125) from 
the APE group using HQET, and M c = 1.33(8) 
from the FNAL group, using an action which in- 
terpolates between relativistic and nonrelativistic 
regions. 
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